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Conjugate gradient (CG) method is well-known for its ability to solve 
unconstrained optimization (UO.) problems. This article presenting a new 
CG method with sufficient descent conditions which improves the former 
method developed by Rvaie, Mustafa, Ismail and Leong (RMIL). The 
efficacy of the proposed method has been demonstrated through simulations 
on the Kijang Emas pricing regression problem. The daily data between 
January 2021 to May 2021 were obtained from Malaysian Ministry of 
Health and Bank Negara Malaysia. The dependent variable for this study 
was the Kijang Emas price, and the independent variables were the 
coronavirus disease (COVID-19) measures (i.e., new cases, R-naught, death 
cases, new recovered). Data collected were analyzed on its correlation and 
coefficient determinant, and the influences of COVID-19 on Kijang Emas 
price was examined through multiple linear regression model. Findings 
revealed that the suggested technique outperformed the existing CG 


algorithms in terms of computing efficiency. 
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1. INTRODUCTION 

The conjugate gradient (CG) is one of the optimization techniques using gradient based methods for 
parameter learning. It is one of the well-known methods to solve unconstrained optimization (UO.) problems. 
The UO. problems usually arise directly in some applications which includes the least-squares problem for 
data-fitting or regression analysis [1]. The structure of the CG method consists of four crucial procedures: 


— Computation of a search direction, dx from x, using (1), 


_ J—8r if k=0, 
E: +prdk-v if k21. (1) 


— Ensure that the direction dẹ denotes a descent direction to guarantee reduction on the objective function 


f. 
— Computation of step-size, a, such that f(x, + a,d,) < f (Xp) is satisfied. The step-size can be computed 
using an exact or inexact line search known as strong Wolfe [2] specified by the rules as in (2) and (3): 
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f(x + a, dx) < fOr) + 0.0001a;,¢7 dx, (2) 
gl (xe + Ody) dy] < —O.1 gi de. (3) 


— Updating the new iterative point, x,4, =X, + kdg. The iteration will be terminated when stopping 
criteria have been satisfied. 

The coefficient 6, in (1) is called the CG coefficient and plays a vital role in the convergence 
behavior. Some of the 6; proposed in the literature were suggested by Hestenes and Stiefel [3], Fletcher [4], 
Polak and Ribiere [5], Polyak [6], Daniel [7], Fletcher [8], Liu and Storey [9], Dai and Yuan [10]. However, 
these classical CG failed to satisfy the sufficient descent condition (S.D.C.) [11] for some inexact line search, 
which guarantee global convergence. The S.D.C. can be represented by (4). 


gkdk < —cllgxll?, k 2 0,c > 0. (4) 


There are various attempts by researchers to rectify this shortcoming, including research by Rivaie, 
et al. [12] who proposed new pp defined as (5). 


T 
RMIL Zk (Sk-8k-1) 
zeree 5 
Pr dt, (dk—1-8k) (5) 


In the same year, the same authors also proposed a similar beta coefficient known as RMIL* in [13]. 


RMIL* _ &4Yk-1 
= ai 6 
By Ida ll? (6) 


Although RMIL/RMIL* has global convergence properties in correspondence to exact line search, 
its numerical output lacks the original method of Hestenes and Stiefel (HS) and Polak-Ribiére-Polyak (PRP). 
Besides, RMIL/RMIL* was less efficient and robust under inexact line search as compared to exact line 
search [14]. Therefore, some variant RMIL/RMIL* have been introduced, for example in [14]-[23]. One of 
the most efficient versions known as RMIL+ method was introduced in research [24] and later analyzed in 
research [25] defined as (7). 


RMIL EEEren) if 0 < |gigr-1] < ligkll? 
Bx t= lldg-4ll? ? = l6k5k-1] = kil» (7) 


0, otherwise. 


However, it was also agreed that the three-term CG algorithms are more efficient, reliable, and 
resilient than traditional CG algorithms [26]. In addition, some three-term CG algorithms may be capable of 
providing reasonable directions that ensure global convergence [27]. The beneficial effects of three-term CG 
have encouraged more investigations of its application to the RMIL method. The following are some of the 
three-term versions of RMIL methods termed as ttRMIL+, ttRMIL and ttRMIL* proposed in [28]-[30], 
respectively. 


ee ifk = 0, 

dawnt = B + pE dk- + fie &k-v ifk2 1. 7 
ee ifk = 0, 

qtm — + pRMILg — qiyo ifk 21. (9) 
a if k=0, 

qttRMILy — P + BRMIbeg | — ae Yk-v İf k21. ii 


In CG methods, where the focus is often limited to their successive iteration formulation, the 
usefulness of the initial direction is yet to be proved. The review paper by [11] pointed out that it is necessary 
to take the direction of negative gradient as the initial search direction. Saha and Nath [31] proved that 
convergence may also be assured by adding parameter 0 < w < 1 to the negative gradient in the following 
manner: 
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_ (“Er + 8x, for k=0 
di T E + Prdk-1» for k> (1) 


Since w is a scalar, then the direction of dọ still follows the direction of the negative gradient. Thus, this 
research proposed further improvement on RMIL method [12] by adopting the three-term parameter and a 
scaled initial direction. The complete definition of the proposed method is: 


8k + Wx, ifk = 0 and w E (0,1), 


RMIL gkdk-1 . 
—g,+ Be dk. +>  Ą=— ifk > 1. 


(12) 
dt, (dea) OF“ 


qNEWRMIL = | 
k = 


To be classified as efficient, a CG technique must at the very least satisfy the S.D.C. and have convergence 
properties. The following theorem shows that our proposed CG method named as NEWRMIL, satisfies the 
S.D.C. NEWRMIL referred to the modified RMIL method. 


Theorem 1.1 (Sufficient descent under exact line search) Assume that a CG method with search direction is 
generated by the NEWRMIL method for all k = 0 under exact line search. Then the proposed method 
satisfies the S.D.C. 


gid, <—cllg|l?, k20c>0. 


Proof. If k = 0, then the initial direction takes dg = —gọ + wgo = —(1 — w)gZo where w E (0,1). Next, 
multiply with g}, 


godo = -(1 — w)gigo 


13 
-(1 - ø)ligoll?. eo 


Since (1 — w) > 0, then the S.D.C. is satisfied for k = 0. The proof is ensued for k = 1 by induction 
T 
_ 8kdk-1 


approach where dp = —g, + BRM“d,_, dg? &k-1: Next, multiply with gf, 
T T RMIL gkdk-1 
8kdk = -8k (2x +k a — atest 8-1) 
= ~Ilgell? — BEM ed ys + gt (5 SS — 2-1) st 
For exact line search, gfd,_, = 0. Thus, gd, = —||g,||?, which implies d, is sufficient descent direction. 


Hence, gd, < —cllg,||? holds true. The proof is complete. 

The remaining parts of this article are arranged as follow: The case problems for this research are 
discussed in Section 2, and the research method are presented in Section 3. Some numerical experiments 
have been done to demonstrate the efficacy of the new method, and results are discussed in 
Section 4. Finally, the conclusion is summarized in Section 5. 


2. DESCRIPTION OF PROBLEM 

The 25" of January 2020 marked the first case of coronavirus disease (COVID-19) in Malaysia 
involving three Chinese citizens and the first local case was detected on the 4" of February 2020. On 30% of 
January 2020, COVID-19 virus was declared as public health concern by World Health Organization 
(WHO). On 11" of March 2020, it was acknowledged as a worldwide pandemic [32]. The spread of 
COVID-19 cases accelerated globally, with high risk of human fatalities, infections, and health issues. The 
rate of disease spread was measured by the index of R-naught, Rọ. Governments, publics and mass media all 
over the world frequently utilize Rọ to raise the alarm, and the Ry higher than 1 is a significant source of 
worry [33]. 

This highly contagious pandemic also harms the global economy [34]. The worldwide uncertainties 
regarding the COVID-19 outbreaks have affected the market dynamics of major commodities, including 
gold. For instance, prior to the current pandemic crisis, the highest gold spot price was USD 60.9807074 or 
RM181.36 per gram on 5“ of September 2011 [35]. Meanwhile, throughout the pandemic, the highest gold 
spot price hit USD 66.286 per gram (RM 278.24) on 7 of August 2020. The price continued to fluctuate, 
reaching a drop to USD 52.233 per gramme (RM 227.61) on 30" of March 2020. 
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There were many studies and discussions on determining the asset property of gold. Baur and Lucey 
[36] investigated the role of gold among investors based on the daily price data of stock markets and closing 
spot gold from the 30" of November 1995 to the 30" of November 2005. They provided evidences that gold 
has a short haven time, which was estimated to be around 15 trading days. They concluded that the 
safe-haven property of gold was only valid after extreme adverse stock market shocks bonds but not during 
the rising of stock markets. Salisu et al. [37] also confirmed the hedging effectiveness of gold with respect to 
risks associated with crude oil, notably during the pandemic era, since gold is a hedge towards geopolitical 
and economic problems. Gold investment during the COVID-19 pandemic appears to be a wise financial 
option, as it will shield investors from market fluctuations in stock and oil markets [38]. The gold market also 
serves as a haven and performs better than U.S. stocks during the pandemic [39]. The empirical study by 
Atri et al. [34] employing the autoregressive distributed lag (ARDL) approach from 23" of January 2020 to 
the 23" of June 2020 concluded that COVID-19 current media coverage, deaths, and infections positively 
affect the gold price. 

The outbreak of the pandemic also increased the volatility of gold return. The COVID-19 impact on 
daily gold spot price from January 2020 to May 2020 was proven in the study carried out by Yousef and 
Shehadeh [40]. The researchers utilized a generalized autoregressive conditional heteroskedasticity 
(GARCH) model with a single dependent variable, the US dollars’ daily gold spot price. The logarithms of 
the cumulative total number of daily cases and the number of fresh possibilities were the independent 
variables. Jianyi et al. [41] investigated the effect of COVID-19 on the gold price by developing a multiple 
linear regression model. They used the Python third-party machine learning algorithm as a framework for 
multivariate linear fitting. The data covers the period from 22™ of January 2020 to 23" of September 2020, 
consisting of the COVID-19 measures (number of confirmed cases, deaths, and cured). The value of R square 
obtained was roughly around 0.926, indicating that the model fits the data well. Motivated by the strong 
potential of gold in protecting hedging properties and its ability in providing a strong financial safeguard 
during any crisis, Malaysian State Government has produce various gold coins to encourage interests 
towards gold investments among the people [42]. Malaysia was the 12 country worldwide to produce its 
first gold coin, the Kijang Emas (KE), introduced by Malaysia’s then-Prime Minister, Tun Dr Mahathir 
Mohamad, on 17" July 2001. 

The goal of this study is to apply the CG algorithm to the KE regression model pricing. The gold 
price model provided by Jianyi et al. [41] was adopted in our KE model, with an extra R-naught reproductive 
index number as the independent variable. As the results, a data set was generated with four independent 
variables linked to COVID-19: new cases (NC), R-naught (Ro), death cases (ND), and new recovered (NR). 
The suggested CG and several current CG techniques were used to generate the regression models. The 
performance of each CG algorithm was compared using error analysis, which includes mean squared error 
(MSE), root mean square error (RMSE), mean absolute error (MAE), and R squared. 


3. RESEARCH METHOD 

The method used in this study is visualized in Figure 1, consisting of three main stages: i) data 
collection, ii) correlation analysis and calculation of coefficient of determination to investigate relationship 
among variables, and iii) multiple linear regression model produced by applying the CG method together 
with the least square approach. The data collection period was done from January to May 2021 and includes 
all independent variables obtained from the ministry of health’s (MOH) daily data. It should be noted that 
these are not the only factors affecting KE pricing, but for this current research, the considered variables were 
limited to the factors related to COVID-19 circumferences. This study’s KE gold bullion price data was 
derived from Bank Negara Malaysia, and it includes the selling price of one gram of KE. To maintain the 
integrity of the five-time series, the discordant sequences caused by weekends or holidays were removed, 
hence yielding a total of 470 sets of data. Table 1 outlines the database’s descriptive analysis of each data 
sets’ minimum, maximum, and mean values. The data was then transformed using Z-score standardization 
method as defined by (15) [43]. 


X—Mean(Xx) 


Z — score = 2 
Standarddeviation(Xx) 


(5) 


The standardized data helps produce dimensionless regression coefficients, especially when dealing 
with iterative optimization algorithms. Based on the conceptual framework depicted in Figure 2, the 
dependent variable KE is modeled through a linear combination of the four independent variables. The 
approximated KE is given by the following multiple linear regression (MLR) model: 


KE; = ao + a1 (Ro)i + a2NCi + a3NR; + a4ND; + Ei. (16) 
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Figure 1. Research methodology 


Table 1. Data description 


Variable Minimum Maximum Mean 
KE (RM per gram) 234.69 268.01 250.5 
4 
Ro 0.81 1.2 1.03 
NC 941 8290 2891 
NR 994 5718 2513 
ND 1 61 14 


Next was the determination of regression parameters ao, 41, A2, A3, A4 through least squares method involving 
the sum of the squares of residuals, E;: 


S, = Xi E? = fea (KE; — ao — a, (Ro)i — @2NC; — a3NR; — a,ND;)°. (17) 


Furthermore, the least-squares approach could be transformed into unconstrained optimization (UO.) 
problem with S, as the objective function: 


min Sp: (18) 


(a0,41,42,43,44) 
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The CG method is prominent in solving UO. problems and the flowchart explaining the procedures 


is presented in Figure 1. The CG method is one of the nonlinear optimization routines starting with a guess 


(0) 0) GO aw 
3 


for the parameters Xp = {Ay ai °,Q;°, j a>. The complete definition of gą corresponding to S, for this 


paper are: 
Srao = —2Y (KE; — ao — a, (Ro)i — a2NC; — a3NR; — a,NDj) 
Sra, = —2Y (Ro) (KE; — ao — a1 (Ro)i — a2NC; — a3NR; — a,ND;) 
Sd, =—2)) NC,(KE; — ao — a,(Ro); — a2NC; — azNR; — a,ND;) (19) 
S,a3  =—2)) NR,(KE; — ao — a, (Ro); — A2NC; — a3NR; — a,ND;) 
Sra, = —2 X ND;(KE; vpe ao aa a,(Ro)i = a,NC; bias a3NR; = a,ND;). 

Consequently, the optimal value of the regression parameters is obtained when ||g;|| < 1076 where |]. || is 


the Euclidean norm and the value 10~° corresponds to the tolerance of accuracy. 


| R-Naught, RO | 


| New cases, NC | 


Kijang Emas price, KE 


| New recovered, NR | 


| New death, ND | 


Figure 2. Conceptual framework 


All CG algorithms were evaluated on the same machine using a strong Wolfe line search to provide 
unbiased numerical results. The test problems that cannot be solved using any methods were labeled as 
°’ FAIL,’ which indicates that the NOI and central processing unit (CPU) time values were unavailable. The 
process was usually halted due to two conditions: it takes a longer time to discover a positive value of step- 
size when iteration is more than 1000 or the method cannot be computed by strong Wolfe. The iteration was 
stopped when it reached the accuracy up to six (6) decimal places, which implies that iterates were already 
adequately close to the solution points. 


4. RESULT AND DISCUSSION 

Simulations were executed in this section to evaluate the performance of the proposed CG technique 
by comparing it with existing CG methods and also the SPSS software using the methodology as described in 
Section 3. The multiple regression model was obtained by employing our proposed three-term CG method, 
named NEWRMIL, compared to RMIL+, ttRMIL+, ttRMIL, and ttRMIL* methods. The parameter value for 
NEWRMIL is set to w = 0.3. These algorithms were carried out using MATLAB R2020b subroutine 
program on AMD Ryzen 5 3550H, Radeon Vega Mobile Gfx, 2,100 Mhz, 4 Core(s), 8 logical processor(s). 

Correlation analysis describes the relations of direct influences between the variables. This paper 
uses the Spearman’s Rank correlation coefficient, —1 < r < 1 [44] to dictate the strength among variables 
such that : 


< 0.35 means weak correlations; 
r = j= 0.68 means strong correlations; (20) 
otherwise means moderate correlations. 


In addition, r = 1 means a perfect relation, and r = 0 means no relation. 
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The result in Table 2 shows the relationship between KE e NC (0.756771) and NC e ND 
(0.88064) are strong positive relationship. KE @ Ry (0.6590331), KE e ND (0.653945), Rọ @ NC 
(0.558109), NR e ND (0.562716), and NR @ NC (0.594737) exhibit moderate positive relationship. 
In contrast, Rọ e ND (0.389516) shows weak positive relationship. To conclude, all relationships are 
significant (p-value < 0.05). 

Apart from that, R? was analyzed to evaluate how much variation the independent variables explain 
the dependent variable, KE. Table 2 shows that the variation in KE is explained by Rg (43.4%), NC 
(57.3%), NR (9.4%), and ND (42.8%), respectively. The variation in Rp is explained by ND (15.2%), and NC 
(31.1%), respectively. The variation in NC is explained by ND (77.6%), whereas the variation in NR is 
explained by ND (31.7%) and NC (35.4%). 


Table 2. Correlation and coefficient of determination analysis 


Relationship _ Spearman’s Rank Correlation Coefficient (r) Significant (p value) Strength R? 

KE © Ro 0.6590331 5.155E-13 moderate 0.4343246 
KE © NC 0.756771 1.12E-18 strong 0.572702 
KE e NR 0.307447 0.002576 weak 0.094524 
KE @ ND 0.653945 8.91E-13 moderate 0.427644 
Ry © ND 0.389516 0.000104 weak 0.151723 
Ro © NC 0.558109 5.09E-09 moderate 0.311486 
NC e ND 0.88064 1.34E-31 strong 0.775527 
NR e ND 0.562716 3.57E-09 moderate 0.31665 
NR e NC 0.594737 2.61E-10 moderate 0.353712 


As shown in Table 2, all the relationships are significant (p-value < 0.05). Thus, the regression 
analysis was carried on to measure the direct effect of the independent variables upon the dependent variable. 
The complete multiple regression output from SPSS is depicted in Table 3. Based on the significant values 
(p-value< 0.05), only Rọ and NC were considered to have impact on KE. 


Table 3. Regression analysis using SPSS 


Variable Unstandardized parameter Standard error Standardized parameter t Sig. 
Constant -1.72E-11 062 000 1.000 
Ro 349 .095 .349 3.685 .000 
NC .480 .178 .480 2.705 .008 
NR -.019 .095 -.019 -.203 .840 
ND .016 .136 .106 -780 437 


Therefore, another linear model was developed using MLR model of Kijang Emas by dropping the NR and 
ND to (16). The complete equation is shown in (21). 


KE; = ao + a, (Ro); + a2 NC;. (21) 


The results of linear regression to demonstrate the benefits of the proposed algorithm is displayed in 
Table 4. By using ten different initial points, the tests which achieved the solution point were considered 
success. Therefore, the outcome of this research will be concluded based on the success percentage of 
problem solved. Results in Table 4 shows that ttRMIL+ failed to solve one problem, but ttRMIL, ttRMIL*, 
and RMIL+ was unable to solve three, four, and ten problems, respectively. However, NEWRMIL has been 
able to solve almost all problems 

The detailed comparison of the efficiency and robustness of all methods was made using the 
performance profile method of Dolan and Moré [45] which is a practical methodology for standardizing the 
comparison of algorithms. The performance criteria used in this research focused on the speed of 
convergence, measured in the number of iterations (NOI), CPU time, and accuracy. Figure 3 shows the 
performance profile of this study based on the output displayed in Table 4. Figure 3(a) and 3(b) specifically 
summarize the plots of performance profile due to NOI and CPU time, respectively. The proportion P of test 
problems for which each method is within a factor t of the best are also shown. The proportion of test 
problems for which a method is the fastest are shown on the left side of the graph. At the same time, the 
percentage of test problems successfully solved by each method, which is also an indicator of their 
effectiveness, are shown on the right side. The curved shape of NEWRMIL is the most top compared to other 
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methods, suggested its effectivity in outperforming all the CG methods in terms of speed and degree of 
robustness. 

The model’s performance was then evaluated in terms of mean squared error (MSE), root mean 
square error (RMSE), mean absolute error (MAE) and R squared (R°). The performance of the method was 
observed through the one with the least MSE, RMSE, MEA with the R? close to 1. Interestingly, the value of 
MSE, RMSE, MAE and R? for all CG methods and the SPSS were MSE ~ 0.342262873, RMSE ~ 
0.585032369, and MAE ~ 0.426035347. All models produced R? ~ 0.654056881, which translated to 
65% of KE explained by the COVID-19 measures (Ro, NC). 


Table 4. Regression coefficients approximates for (21) 


Initial point Ag ay az CPU time NOI 
SPSS NIL -1.64E-11 0.343743 0.564925 
NEWRMIL (5,1,0) -1.4916E-09 0.343743421 0.564924647 0.0003495 9 
RMIL+ FAIL 
ttRMIL+ -1.14937E-10 0.343743429 0.564924643 0.0027048 15 
ttRMIL FAIL 
ttRMIL* FAIL 
NEWRMIL (4, 3, -2) -3.67718E-10 0.343743421 0.564924647 0.0003403 9 
RMIL+ -2.25769E-10 0.343743426 0.564924645 0.0003854 10 
ttRMIL+ -2.251E-10 0.343743426 0.564924645 0.0004039 10 
ttRMIL -1.9463 1E-09 0.343743424 0.564924645 0.000364 
ttRMIL* -1.94632E-09 0.343743424 0.564924645 0.0003956 9 
NEWRMIL (3, 3, -1) -5.25E-10 0.34374342 0.564924648 0.0003863 11 
RMIL+ 4.62E-09 0.343743422 0.564924646 0.0005053 12 
ttRMIL+ -2.67E-11 0.343743422 0.564924647 0.0014332 16 
ttRMIL -1.62E-09 0.343743426 0.564924644 0.0005956 14 
ttRMIL* -1.62E-09 0.343743426 0.564924644 0.0004407 14 
NEWRMIL (2, 2, -1) -5.57267E-10 0.343743419 0.564924649 0.0003861 11 
RMIL+ FAIL 
ttRMIL+ -1.30976E-09 0.343743424 0.564924647 0.0020928 18 
ttRMIL 2.50173E-10 0.343743426 0.564924647 0.0003834 11 
ttRMIL* 2.51673E-10 0.343743426 0.564924647 0.0003621 11 
NEWRMIL d, 1,5) 8.67E-10 0.343743417 0.56492465 1 0.0004388 13 
RMIL+ 7.67E-11 0.343743421 0.564924648 0.000769 14 
ttRMIL+ FAIL 
ttRMIL 1.57E-09 0.34374342 0.56492465 0.0004625 13 
ttRMIL* 1.57E-09 0.34374342 0.56492465 0.0005417 13 
NEWRMIL  (-1,-1, 4) 8.37E-10 0.343743421 0.564924647 0.0004226 12 
RMIL+ -2.32E-11 0.343743421 0.564924647 0.0003932 12 
ttRMIL+ -1.90E-12 0.343743421 0.564924648 0.0020654 18 
ttRMIL 5.90E-11 0.343743419 0.564924647 0.000479 15 
ttRMIL* -1.25E-09 0.343743417 0.564924649 0.0007025 16 
NEWRMIL (-2, 0, 4) -2.48611E-09 0.343743415 0.564924653 0.0009874 18 
RMIL+ FAIL 
ttRMIL+ -1.77212E-10 0.343743416 0.56492465 0.0004706 13 
ttRMIL FAIL 
ttRMIL* FAIL 
NEWRMIL (-3, -2, 3) -1.67511E-11 0.343743422 0.564924647 0.0003677 10 
RMIL+ -2.58514E-12 0.343743416 0.56492465 0.0005915 12 
ttRMIL+ 2.0645E-09 0.343743419 0.564924648 0.0004075 11 
ttRMIL FAIL 
ttRMIL* FAIL 
NEWRMIL (-4, -1, 3) 6.10535E-10 0.343743422 0.564924646 0.0004766 11 
RMIL+ 1.96786E-12 0.343743421 0.564924647 0.00042 12 
ttRMIL+ 2.68306E-12 0.343743421 0.564924647 0.0004491 12 
ttRMIL 6.5898E-10 0.34374342 0.564924649 0.0006077 14 
ttRMIL* 6.59261E-10 0.34374342 0.564924649 0.0004684 14 
NEWRMIL (-5,-1, 1) -2.94735E-09 0.343743417 0.56492465 0.0004612 12 
RMIL+ FAIL 
ttRMIL+ 3.36917E-10 0.343743415 0.564924651 0.0058446 33 
ttRMIL 2.33855E-10 0.34374342 0.564924648 0.0005962 11 
ttRMIL* FAIL 


The overall F statistic is 86.02451223 and the p — value = 1.05856E — 21 < 0.05. Meanwhile, 
t-statistics for KE, Rọ and NC are significant at the 0.05 level. Based on this justification, only Rọ and NC 
impacts KE. This indicates that as the rate of Rọ and the number of new COVID-19 cases rises, the price of 
KE will increase as well. The related regression output is depicted in Table 5. 
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Figure 3. Performance profile based on (a) NOI, and (b) CPU time 


Table 5. Regression analysis using SPSS for two independent variables 


Coefficients Standard Error t Sig. 
Constant -4.18748E-15 0.061328056 -6.828E-14 1 
R-naught 0.343743421 0.07430623 4.626037682 1.22988E-05 
new cases 0.564924647 0.07430623 7.602655183 2.52577E-11 


5. CONCLUSION 

The primary emphasis of this paper is to present a new CG algorithm, NEWRMIL, and its 
applicability in solving linear regression models. The current pandemic crisis directed the author’s interests 
to study the direct impacts of COVID-19 on the price of Kijang Emas (KE) through a multiple linear 
regression model. Estimation of parameters in linear regression has been carried out using the least- squares 
method and CG algorithm. The NEWRMIL method was compared with a set of well-known variants of 
RMIL: RMIL+4, ttRMIL+, ttRMIL, and ttRMIL* methods. Numerical analysis proves that the NEWRMIL 
method has the highest NOI and CPU time efficiency among the CG methods used in this study. Therefore, 
the NEWRMIL method is deemed to be efficient, theoretically and practically. The regression analysis 
results revealed that RO and the number of COVID-19 daily cases significantly affected the price of Kijang 
Emas. Furthermore, the results from the proposed CG method are comparable to SPSS with good 
performance of MSE, MAE, AND R2. In a nutshell, NEWRMIL method is considered as reliable, with 
regard to its practicality, useful estimates and reasonable accuracies. 
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